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Abstract 

The transverse Nernst Ettingshausen (N — E) effect and electron mobility in Pbi_ 2 ;Sn 3 ;Se alloys 
are studied experimentally and theoretically as functions of temperature and chemical composition 
in the vicinity of vanishing energy gap Eg. The study is motivated by the recent discovery that, 
by lowering the temperature, one can change the band ordering from trivial to nontrivial one 
in which the topological crystalline insulator states appear at the surface. Our work presents 
several new aspects. It is shown experimentally and theoretically that the bulk N — E effect 
has a maximum when the energy gap Eg of the mixed crystal goes through zero value. This 
result contradicts the claim made in the literature that the N — E effect changes sign when the 
gap vanishes. We successfully describe dc transport effects in the situation of extreme band’s 
nonparabolicity which, to the best of our knowledge, has never been tried before. A situation is 
reached in which both two-dimensional bands (topological surface states) and three-dimensional 
bands are linear in electron k vector. Various scattering modes and their contribution to transport 
phenomena in Pbi-ajSn^iSe are analyzed. As the energy gap goes through zero, some transport 
integrals have a singular (nonphysical) behaviour and we demonstrate how to deal with this problem 
by introducing damping. 
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I. INTRODUCTION 


Narrow gap semiconductors have been for many years subject of intense experimental 
and theoretical studies in view of their interesting properties and important applications 
In recent years they have become once again a source of excitement due to the dis¬ 
covery of topological insulators [ 1 ] . The topological boundary states have been observed in 
bulk compounds 612803, 612X03, bulk alloys 6ii_3.Sb3;[5], two-dimensional quantum wells of 
HgTe/Hgi_3,Cda;Te | 1 ] and, most recently, in bulk Pbi_3;Sna;Se alloys |6], bulk SnTe [ 7 j and 
bulk Pbi-xSui^Te alloys |H]. The latter were called topological crystalline insulators (TCI), 
because, in contrast to canonical topological insulators, specihc crystalline symmetries war¬ 
rant the topological protection of their metallic surface states [ 9 IITT]. Since the IV-VI lead 
chalcogenides are characterized by strong temperature dependence of their band structures, 
it is possible to reach proper band ordering for a suitably chosen chemical composition x by 
varying the temperature. For the Pbi_a;Sna;Se system at high temperatures the band order¬ 
ing is called trivial (Lg band above Lq band) and no TCI state occurs. As the temperature 
is lowered for the properly chosen x, one can reach the vanishing band gap and then arrive 
at the inverted nontrivial ordering band above Lq band) in which TCI state can occur. 
Such transition was demonstrated with the use of angle-resolved photoemission spectroscopy 
in Ref. [B] . Generally speaking, the character of topological insulators is manifested in both 
their surface and bulk properties. 

It has been recognized from the sixties that in the alloys Pbi-^Sn^^Se one can reach van¬ 
ishing energy gap by changing the temperature. An early demonstration of this possibility 
was provided by Strauss |T2] who used for this purpose the optical transmission measuring 
the gap on both sides of the band ordering. Also laser emission proved useful in this re¬ 
spect [ 13 ]. It was shown that the electric resistivity and the Hall coefficient depend on the 
band ordering and can be used to obtain information on the transition temperature even 
if the Fermi energy is quite high in either the valence or conduction band na. Quite re¬ 
cently, optical response with appropriate analysis was used to determine the band inversion 
temperature for Pbo,77Sno.23Se [TSlITT] . 

In recent years the thermomagnetic Nernst-Ettingshausen effect experiences a real revival 
in investigations of semiconductors [131 [HI, graphene [Ml El] and high-Tc superconductors 
[221E 3 ]. In the present work we undertake a bulk transport study of PbSe and Pbi_a;Sn3;Se 


system for 0.25 < x < 0.39. The above ternary alloys can reach the vanishing band gap 
in the available temperature range. We investigate experimentally and theoretically the 
conduction electron density, mobility and the N — E effect. Since our study concentrates 
on small gap values, particular features of the band structure are crucial and it turns out 
that the description of transport effects for the vanishing gap poses nontrivial theoretical 
problems. In fact, although we describe bulk properties, as the gap goes to zero we deal 
with linear energy bands which is a completely new situation for the transport theory. We 
demonstrate for four Pbi_ 3 ;Sn 3 ,Se samples that, as the gap goes to zero, the N — E effect 
reaches a maximum. This contradicts the claims made in the literature, see Ref. [18]. All 
in all, our results conhrm the conclusions of Ref. [6] concerning the transition from trivial 
to nontrivial band ordering as a function of temperature. Our analysis of scattering modes 
and their relative importance in Pbi_a;Sna,Se near the trivial-nontrivial transition of the 
band ordering will help further investigations of topological crystalline insulators and other 
topological materials. 

II. EXPERIMENT 

Single crystals of Pbi_a,Sna,Se (0< x < 0.39) were grown by self-selecting vapour growth 
technique [SUES]. Owing to the peculiarity of the method, where near equilibrium thermo¬ 
dynamic conditions are kept, we have obtained high quality compositionally uniform large 
monocrystals with natural (001) facets (typical dimensions: 1x1x1 cm^). Crystal com¬ 
positions were determined by Energy Dispersive X-ray Spectroscopy offered by Scanning 
Electron Microscope Hitachi SU-70. Molar fractions were taken as averages from the scan 
covering lx 1.5 mm^ of the sample surface. The accuracy of determination of the chemical 
composition is better than 0.005 molar fraction. 

Samples for measurements were cleaved with razor blade along (001) planes in the form of 
rectangular parallelepipeds, with dimensions 1.5x3x10 mm^. All samples exhibited metal¬ 
lic behaviour of resistivity with almost temperature independent carrier concentration. The 
Hall effect measured electron density in PbSe to be 7T0^® cm“^ and the mobility 35000 
cm^/Vs at liquid helium temperature. The corresponding values for ternary compounds 
Pbi_a;Sna;Se are included in Table 1. For thermoelectric measurements, the samples were 
thermally anchored with silver epoxy to copper cold finger of continuous flow helium cryo- 


stat. A small SMD resistor was glued to the free end of the sample for use as a heater. 
Copper potential leads were attached across the sample with silver paint. Temperature and 
temperature gradient were measured with two calibrated subminiature GaAlAs diodes by 
Lake Shore glued to the sample. For all temperatures, the temperature gradient was kept 
within 2-5 percent of the average temperature. An external magnetic held at given temper¬ 
ature, with proper temperature gradient steady developed, was swept from -0.5 to 0.5 Tesla 
using standard resistive magnet. The magnetic held was directed perpendicularly to the 
sample and the temperature gradient. In this conhguration the transverse N — E electric 
held develops across the sample: 


-En-e = -Pn-e ■ VT X B , (1) 

where Pn-e is the N — E coefficient, T is the temperature and B is the magnetic held. 
The transverse N — E ehect is a thermoelectric analogue of the Hall ehect. The potential 
diherence Vn-e=w-En_e (w is sample’s width) was measured by Keithley nanovoltmeter 
2182a. The potential was linearly dependent on the magnetic held for all applied helds. 
Main uncertainties in the absolute Pn-e data are systematic and arise mainly from hnite 
size of the thermometers contacts to the sample. We estimate their contribution to the 
uncertainty of the temperature gradient to be less than 15 percent. All other possible errors 
are small in comparison. 


III. THEORY 

The conduction band of PbSe and Pbi_j,Sna:Se alloys consists of four ellipsoids of revo¬ 
lution with minima at the four L points of the Brillouin zone. As a consequence, one deals 
with the longitudinal m*i and transverse ehective masses, the anisotropy in PbSe at 0 
K is m'^/m^= 1.7. In the description of electron scattering one deals with the density-of- 
states mass m*^ = and in the description of mobility with the conductivity mass 

m* = + 2ml~^). In PbSe and Pbi_a;Snj;Se these masses diher little from each other, 

so we use an approximate spherical value of tUq = 0.048mo at 0 K at the band edge (mo 
is free electron’s mass). This means that we use the standard dispersion relation for the 



two-band Kane model (the zero of energy is chosen at the band edge) 
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This gives the energy-dependent mass, see [23 

m*(T) = m* (^1+ 2^^^ , (3) 

It is seen from Eqs. (2) and (3) that, as the gap vanishes, we deal with the linear energy 
dispersion E(k) i.e. with, so called, massless Dirac fermions. In principle the k-p theory is 
valid at T = 0, bnt here it is known to work also for higher temperatnres. Using experimental 
data on the composition dependence of the band gap of Pbi_ 3 ,Sna,Se at varions temperatnres, 
see Refs. [IlE2l|27|, we establish the following phenomenological dependence of Eg (meV) 
on X and T 

Eg{T, x) = 126 - 802 ■ a; + [400 + 0.258 ■ , (4) 


In the two-band k-p model the effective mass is proportional to the gap and inversely 
proportional to the matrix element of momentnm sqnared. The temperatnre dependence of 
the mass at the band edge of Pbi_a;Sna;Se is 
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mo 3.0546 + 1.8626 ■ 10-3 ■ T 


(5) 


In order to consistently describe transport effects we nse the formalism developed for non¬ 
parabolic energy bands, see [26]. According to this scheme, statistical and transport qnan- 
tities are given by the integrals 

< A ’ (6) 

where /o is the Fermi-Dirac distribntion fnnction depending on the Fermi level Ef. Since the 
derivative dfQ{S)/dS does not vanish only in the limited range of energies aronnd the Fermi 
level, the integrals (6) are not difficnlt to compnte. If the electron gas is strongly degenerate, 
there is dfo{S)/dS = —6{S — Ef) and the integrals (6) are eqnal to the integrands taken at 
the Fermi energy. The free electron density in the band is 
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where iV^ = 4 is the number of equivalent ellipsoids. The electron mobility is ^{£) = 
eT{£)/m*{£) , in which r is the relaxation time and m*{£) is given by Eq. (3). The electric 
conductivity is a = eNjl, where the average electron mobility is 
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An experimental determination of free electron density is nontrivial for our samples since it 
involves the Hall scattering factor 
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which for materials with vanishing energy gap should be carefully evaluated, see below. The 
N — E effect, which is the main subject of our interest, is described in Ref. 
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where: ks - Boltzmann constant and For the complete degeneracy of electron 

gas the difference in the parenthesis is zero and the N — E effect vanishes. 

In order to describe electron scattering mechanisms one needs to know the electron wave 
functions for the conduction band. These are taken in the form of true Bloch states with the 
periodic amplitudes depending on the pseudo-momentum Kk and energy, see [28]. Because 
of the narrow energy gap and fairly strong spin-orbit interaction the wave functions are in 
general linear combinations of conduction and valence periodic states as well as mixtures 
of spin-up and spin-down components. The wave functions for the conduction band in the 
spherical approximation are, cf. 




Vl - LbZ + VZ^iR 
k 


t 


' _ u 

\/l — LdX^ -|- vL—iR 
k 


; , 


(11) 


^7 = 


Vl - LhZ + y/Z^iR 
k 


i + 


\/l - LdX_ + VZ—iR 
k 


t, 


( 12 ) 


where L = £/{2£ + Eg) and k± = k^ A: iky, while Z, X± = (X± \Y)/^/2 and R denote the 
periodic amplitudes of Luttinger-Kohn functions taken at the L points of the Brillouin zone. 
The normalization coefficients satisfy the condition IZ + d?=l. Plus and minus signs in the 




















overscripts of the wave functions symbolize the effective spins up and down, respectively, 
while the arrows on the RHS symbolize spin-up and spin-down functions quantized on the 
z direction. The electron energies E are counted from the conduction band edge. 

The scattering probabilities have been calculated for various scattering modes by comput¬ 
ing matrix elements of corresponding perturbing potentials. The total scattering probability 
is a sum of separate probabilities, which amounts to calculating the total relaxation time 
according the well known formula 

i 

where r* describe relaxation times for specihc (independent) modes. Below we enumerate the 
relevant scattering modes, quote the involved material parameters and mention importance 
of various modes for the electron mobility and N — E effect in Pbi_j,Sna,Se for the chemical 
compositions 0.25 < x < 0.39 of our interest. The description includes band’s nonparabol- 
icity in the energies and wave functions, as well as screening of long-range interactions by 
the electron gas. 

1) Polar scattering by optic phonons (OP). This mode is important at higher tempera¬ 
tures. It is determined by the static and high-frequency dielectric constants Eq =234 and 
£oo = 28.7 at T = 0. Their temperature variations can be found in [H 121]. The mode 
is characterized by nonelastic scattering processes which complicates its description in 
relaxation time approximation. This mode contains no adjustable parameters. Taking into 
account the screening of the interaction by the electron gas one obtains for the relaxation 
time, cf. [2H] 
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in which poo = 4/c^Aoo^, where Aoo is the screening length for ^oo- 

2) Nonpolar scattering by optic phonons (NOP). This mode is determined by electron- 
optic-phonon deformation potential interaction. Due to the interband k-p mixing it involves 









both the conduction and valence deformation potential constants, and which are 
treated as adjustable parameters, cf. ISI- We assumed that and are equal. A 
correction due to the nonelasticity is included, see EH. The relaxation time for this mode is 
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where 
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Here hwop is the energy of the optical phonon and Oq is the lattice constant. The correcting 
term T^or = (expz — exp 2 ;), where z = S/ksT. 

3) Scattering by acoustic phonons (AC). It involves conduction and valence acoustic 
deformation potentials E^^ and E^^ treated as adjustable parameters. We assumed that E^^ 
and E^^ are equal - a good approximation for mirror-like conduction band - valence band 
symmetry in IV-VI semiconductors. The relaxation time for this mode is 
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where Vav is the averaged sound velocity and g is the crystal density. 


4) Scattering by ionized defects is due to electrostatic interaction between electrons and 
charged defects in the crystal. In Pbi_a;Sna;Se each native defect furnishes two free electrons 
(Se vacancies) or two free holes (metal vacancies) [1]. In general, the defect potential is of the 
form V=Vc + V^,., where Vc is the Coulomb interaction and V^^ symbolizes the short range 
interaction related to size of the defect. The Coulomb interaction is negligible in lead chalco- 
genides due to very high value of the static dielectric constant Eq. Thus, one is left with the 










short-range contribution to electron scattering. The relaxation time for this mode is, see [32] 
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where Nd is the concentration of ionized defects, while A = < R\Vsr\R > and B=< 
A|14r|-A > are the matrix elements of the short-range potential Vsr for the conduction and 
valence bands, respectively. The elements A and B are treated as adjustable parameters 
and were assumed to be equal. 

5) Alloy disorder scattering (AD), that appears only in the ternary alloys, is due to the 
fact that the Vpb and Vsn atomic potentials are not the same. This results in perturbations 
of crystal periodicity and, consequently, in electron scattering, see [331 El]- The disorder 
scattering is important at high values of x and low temperatures. The relaxation time for 
this mode is, see [33] 
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where U^d matrix element of the potential difference Vpf, — Vsn for the conduction 

band and D is the volume of unit cell. Further 
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where U^d is the matrix element of the potential difference Vp^-sn for the valence band. 
According to the theory, see [33], there exists a relation 
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Knowing the gaps of both materials and htting the value of one automatically obtains 
the value of U^d- 









6) Scattering by charged dislocations (DIS). This mode is due to repulsive interaction of 
electrons with dislocation lines which, forming acceptor centers, attract conduction electrons 
and become negatively charged, see [35]. The mode depends strongly on dislocation density 
(which is an adjustable parameter) and rather weakly on the temperature. The relaxation 
time for this mode is 


+ ie/E,) 

Tdrs 

where / is the fraction of filled traps, a is the lattice constant, Aq is the screening length for 
Eq and Ndis is the dislocation density. 

In Table 1 we quote the employed values of material parameters. After the adjustment 
they are kept the same for all four samples. The density of dislocations Ndis and the damping 
constant G are adjusted for each sample separately. 

Trying to compare experimental transport data with the theory one has to connect the 
measured quantities with the calculated ones. To work out this question we proceed in the 
following way. We begin with the experimentally determined Hall constant Rh for a given 
temperature T from which we determine the Hall electron density Nh = l/eRn- The real 
free-electron density is = ArNn, where Ar is the Hall scattering factor. From the given 
value of Nh we compute the first value of the Fermi energy Ej and use Ef to compute 
the first value of mobility Jt using Eq. (8) with the above mentioned scattering modes, and 
finally the first value of Aj. given by Eq. (9). Employing this value of Aj. we determine the 
first value of N which concludes the first round of iteration. For this value of N we compute 
the second value of Ef, etc, all the way to the second value of N. We terminate the iteration 
when the (n+l)th value of N is practically the same as the nth. This procedure determines 
the final values of N, Ef, A^ and p for the given temperature. These values can then be used 
to calculate the N — E coefficient given by Eq. (10). The adjustable parameters mentioned 
above are chosen to obtain an overall best agreement between the experiment and theory 
for Jl and Pn-e for our four samples at all available temperatures. 

When employing the above scheme we have found that the computed quantity containing 
the integral < /i^ >, that is Ar and Pn-e, (see Eqs. (9) and (10)), have a strong and narrow 
peak at the temperature Tc for which Eg = 0. Such peak has no physical meaning and, as it 
can be seen in Fig. 1 for Pn-e, h is not observed experimentally. Also, there is no reason to 
expect a sharp peak of the scattering factor Ar because it would lead to a sharp peak of the 



free electron density N withont a physical reason. However, since onr Pbi-ajSnj-Se samples 
are certainly not homogeneons, i.e. they have somewhat different chemical compositions x 
at varions parts, one can not expect to have Eg = {) m the whole sample at one temperatnre 
T. We simnlate this nonhomogeneity by nsing the well known mathematical measnre to 
avoid singnlarities in resonances. Thus, we introduce damping by replacing in all formulas 
the value of Eg by Eg + iG and take real values of the resulting expressions. By adjusting the 
damping constant G we can bring both the scattering factor and the N — E coefficient 
into a reasonable and experimentally observed behaviour, see below. 

IV. RESULTS AND DISCUSSION 

Figure 1 shows our main experimental results on the N — E effect in PbSe and four 
Pbi_a;Sn 2 ;Se samples in which the gaps Eg in Pbi_a;Sn 3 ;Se alloys can go through zero as 
functions of temperature. Our absolute experimental values of Pn-e are in general agree¬ 
ment with those measured on other narrow-gap semiconducting materials [m [36]. The 
vertical thin arrows indicate critical temperatures which, according to the dependence 
(4), correspond to = 0 for the indicated chemical compositions. It can be seen that 
the maximum of Pn-e for each sample corresponds to the critical temperature Tc. In other 
words, a maximum of Pn-e indicates that the gap goes through zero value. On the other 
hand, in PbSe there is no possibility of reaching vanishing gap and one observes no maxi¬ 
mum. The observed values of Tc, are compared to the line drawn on the basis of Eg{T,x) 
given in Eq. (4) and presented in Fig. 2. The agreement between the two is very good. 

This result contradicts the claim made in Ref. na that the temperature at which the 
band gap vanishes corresponds to the change of sign of Pn-e- In particular, in the samples 
with X = 0.25 and x = 0.277 the coefficient Pn-e does not change sign at all but the gap 
goes through zero and the maxima are well observed. In principle, the N — E effect should 
go to zero as the temperature goes to zero, but this is experimentally difficult to achieve 
because it requires a very small temperature gradient. 

In order to understand and appreciate contributions and relative importance of various 
scattering modes in Pbo. 75 Sno. 25 Se, which are used in the description of N — E effect, we show 
in Fig. 3 an example of experimental mobility Jl for this sample, compared with the theory 
for the total mobility and partial mobilities related to single modes. The values of adjusted 


material parameters are given in Table 1. As indicated above, there is JI = n{H)/Ar, so that 
we have to go here through the calculation of Ar as well. As the temperature goes from 0 
to 300 K the calculated value of Ar goes smoothly from 0.85 to 1.55. It can be seen that, 
with the indicated values of material parameters for various scattering modes, the overall 
description of mobility is quite satisfactory. It follows from the hgure that the dominant 
scattering mechanism is due to alloy disorder. At higher temperatures the polar and nonpolar 
optical, as well as acoustical phonon modes become important. For the assumed low density 
of linear dislocations Ndis = 10 ® cm“^ the corresponding partial mobility is too high to be 
seen in the hgure. 

In Fig. 4 we plot theoretical behaviour of Pn-e for Pbi_a;Sna;Se having four different 
chemical compositions x. All calculations are carried for small density of dislocations Ndis = 
1 ■ 10® cm“^, no damping {G = 0) and hxed values of other material parameters. In addition, 
the electron density A^ = 1 ■ 10^® cm“® is kept the same for all samples. This somewhat 
hypothetical calculation is instructive for several reasons. It shows that for a; = 0 and x 
= 0.15, when no = 0 occurs, the behaviour of Pn-e is very hat. On the other hand, 
for higher values of x, where the gap vanishes as a function of T, the N — E ehect has a 
maximum at the critical temperature T^.. It is seen that, when no damping is introduced, 
theoretical Pn-e coefficients have high and very sharp peaks not observed experimentally. 
This, as we explained previously, is a result of the ’’explosive” behaviour of the < > 

integral at = 0 , which appears in formula ( 10 ) for Pn-e- 

Figure 5 shows two calculations of the N — E coefficient for the Pbo. 75 Sno. 25 Se sample, 
carried out with the same values of adjustable material parameters. The dashed line indicates 
a calculation without damping, i.e. with G = 0. It gives a sharp peak at the temperature 
Tc = 141 K at which Eg = 0. This peak is related to the integral jj? discussed above and it 
is not observed experimentally. For this reason we use the damping procedure mentioned in 
the preceding section with the adjusted value of G = 20 meV. This gives the result indicated 
by the solid line. It is seen that now there is no sharp peak and the theory describes quite 
well the experiment shown in Fig. 1 with the exception of low temperatures. Finally, we 
indicate a marked difference of theoretical results for Pn-e calculated for the sample x = 
0.25 with the use of G = 0, as seen in Figs 4 and 5. This difference is due to largely 
different dislocation densities assumed for both calculations: Ndis = 10® cm“^ in Fig. 4 and 
Ndis = 10® cm“® in Fig. 5. 



Finally, Fig. 6 shows theoretical dependence of Pn-e versus temperature, calculated for 
all four Pbi_a;Sna;Se samples of our interest. In the description we use for all samples the same 
values of material parameters, excepting the density of defects, the density of dislocations, 
and the damping constant G. The latter are adjusted for each sample and their values are 
used also in the description of average mobilities, as illustrated for the sample with x = 0.25 
in Fig. 3. In the same hgure we reproduce for comparison the experimental data shown in 
Fig. 1 with the same vertical arrows indicating the critical temperatures at which Eg = 0. 
The vertical arrows drawn at maxima of Pn-e denote the critical temperature Tc at which 
Eg=0. The same arrows are transferred to Fig.l for comparison with experimental data. 
It is seen that the calculated and experimental maxima occur at the same temperatures 
within several percent accuracy. The general experimental and theoretical temperature 
behaviour near the maxima as well as the absolute values of Pn-e agree quite with each 
other. By comparing Figs. 4 and 5 it becomes clear that the maximum oi N — E effect 
is directly related to the gap going through zero. This is a general feature of the theory 
which should find applications for other zero-gap semiconductor and semimetal systems. On 
the other hand, the low temperature behaviour of the samples with x=0.325 and x=0.39 
remains unexplained according to the presented theory. It should be observed that the 
unexplained region corresponds to the nontrivial band ordering in which the topologically 
protected surface states appear, so we deal with a two-channel transport. According to the 
formulas for two-channel mixed conductivity, the N — E coefficient depends on weighted 
contributions of both conducting channels and difference of their thermoelectric powers EH. 
Consequently, it is possible that the low temperature behaviour is modihed by the topological 
surface states which are not accounted for in our theory. It should be mentioned that higher 
values of the damping constant G result in wider and lower theoretical maxima of Pn-e. 
All in all, we achieve a good description of the N — E effect and electron mobility for four 
investigated samples in the critical range of small forbidden gaps in which the transition 
between trivial-nontrivial band ordering takes place. Our analysis has general significance 
for zero-gap systems and, in particular, it should be useful for a description of electron 
transport in topological crystalline insulators and in recently discovered three-dimensional 
topological Dirac semimetals - bulk analogues of graphene [Ml EH] ■ 


V. CONCLUSIONS AND SUMMARY 


The main conclusion of our work is that a maximum of the electron Nernst-Ettingshausen 
effect, measured as a function of decreasing temperature in Pbi_a;Sna;Se alloys, indicates a 
critical temperature Tc corresponding to the zero energy gap and transition from trivial to 
nontrivial band ordering in terms of the topological crystalline insulators. The above prop¬ 
erty allows one to determine experimentally critical temperatures by means of transport 
data. In particular, we checked theoretically that the maximum oi N — E effect does not 
depend on the carrier density in the sample. Electron mobility and N — E effect are mea¬ 
sured and successfully described for four Pbi_j;Snj;Se samples with 0.25< x <0.39 in the 
vicinity of = 0, establishing the dominant scattering modes. It is demonstrated that 
the description of Hall scattering factor and N — E effect presents theoretical problems as 
the gap goes to zero and the bands become linear in the wave vector. The latter have been 
overcome introducing damping procedure. We emphasize that our conclusion concerning 
the maximum oi N — E effect is a result of the detailed theoretical treatment. Our analy¬ 
sis of electron transport phenomena in the proximity of vanishing energy gap should be of 
value for other temperature and composition driven topological crystalline insulators and 
three-dimensional Dirac semimetals. 
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TABLE I; Carrier densities and mobilities, as determined from Hall measurements of Pbi_ 2 ;Sna;Se 
alloys and material parameters used in the calculations. Scattering parameters are adjusted and 
kept the same for all samples. For alloy disorder: -1.25T0“^^ eV cm^, -0.477, 

for short range potential A =1.6T0“^^ eV cm^, for acoustic deformation potential = 35 eV, 
for optic nonpolar deformation potential; = 60 eV. Dislocation density N^is and damping 


constant G are adjusted for each 

sample. 
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FIG. 1: Experimental values of the transverse N — E coefficient, Pn-Ej versus temperature for 
PbSe and four samples of Pbi_xSna;Se with different chemical compositions. Thin vertical arrows 
point to critical temperatures for which the band gap vanishes according to Eq. (4). It is seen 
that Pn-e has a maximum when the gap vanishes. 














FIG. 2: Experimental temperatures of the maxima of Fn-e (Tc) taken from Fig. 1, for four 
investigated samples versus chemical composition x, compared with the line drawn according to 
Eq. (4). 
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FIG. 3: Electron mobility Jl versus temperature for Pbo. 75 Sno, 25 Se sample. Experiment - dots, 
theory - solid line. Theoretical partial mobilities related to separate scattering modes are also 
indicated. Notation: AD - alloy disorder, OP - polar optical phonons, NOP- nonpolar optical 
phonons, AC - acoustic phonons, SR - short range potential. Employed material parameters are 
given in Table 1. 
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FIG. 4: Calculated Pn-e coefficient for PbSe and four Pbi_a;Sn 3 ;Se samples versus temperature 
keeping the same material parameters, the same electron density N and excluding damping. The 
assumed density of dislocations N^ig = 10® cm“^ is very small. For zero damping, the Pn-e 
coefficients have sharp and narrow maxima at the critical temperatures Tc for which the gaps 


vanish. 
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FIG. 5: The N — E coefficient versus temperature calculated for Pbo. 75 Sno. 25 Se sample. Material 
parameters are given in Table 1. The dashed line is computed without damping, i.e. with (7 = 0. 
The solid line includes damping [G = 20 meV). 
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FIG. 6: Theoretical and experimental Pn-e coefficients versus temperature for four investigated 
Pbi_a;Sna;Se samples. Thick solid lines - theory with the use of material parameters given in Table 
1, open circles - experiment. Thin vertical arrows are drawn at the critical temperatures Tc for 
which the corresponding gaps vanish (the same arrows are also drawn in Fig. 1). For temperatures 
near maxima the theory describes very well the experimental data both in terms of shape and 
absolute values. 












